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New exact interior solutions to the Einstein field equations for anisotropic 
spheres are found. We utilise a procedure that necessitates a choice for the 
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\ energy density and the radial pressure. This class contains the constant 

density model of Maharaj and Maartens (Gen. Rel. Grav., 21, 899-905, 
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1989) and the variable density model of Gokhroo and Mehra (Gen. Rel. 
Grav., 26, 75-84, 1994) as special cases. These anisotropic spheres match 
smoothly to the Schwarzschild exterior and gravitational potentials are 
well behaved in the interior. A graphical analysis of the matter variables 
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is performed which points to a physically reasonable matter distribution. 
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1 Introduction 



In recent years a number of authors have studied solutions to the Ein- 
stein field equations corresponding to anisotropic matter where the radial 
component of the pressure differs from the angular component. The grav- 
itational field is taken to be spherically symmetric and static since these 
solutions may be applied to relativistic stars. A number of researchers have 
examined how anisotropic matter affects the critical mass, critical surface 
rcdshift and stability of highly compact bodies. These investigations are 

dtannnn 

contained in the papers la, 13, HfJ, HJ, Ha, Ha, I2H, among others. Some 
researchers have suggested that anisotropy may be important in under- 
standing the gravitational behaviour of boson stars and the role of strange 
matter with densities hi ghe r than neutron stars. Mak and Harko |15| and 



Sharma and Mukherjee 



22( suggest that anisotropy is a crucial ingredient 



in the description of dense stars with strange matter. 



In this paper our objective is to generate a new class of exact solutions to 
the Einstein field equations corresponding to a physically reasonable form 
for the energy density. A particular motive is to find simple analytic forms 
for the gravitational and matter variables so that the physical interpretation 
of the model is simplified. Often the solutions arc presented in terms of 

OQ 

special functions or a numerical approach is required [a, |7j ■ We hope that 
our results in terms of elementary functions will assist in the analysis of 
gravitational behaviour of compact objects, and the study of anisotropy 
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under strong gravitational fields. In section [3 we develop the anisotropic 
stellar model and present the relevant field equations. A particular form for 
the energy density is chosen in section[31 and the Einstein field equations are 
integrated. Special cases of physical interest are isolated from the general 
solution in section 0] Some physical features of the anisotropic star are 
briefly considered in section^ In scction[S]we demonstrate that our model 
yields surface redshifts and masses that correspond to real sources, and 
make a few concluding remarks and suggestions for future research. 



2 The anisotropic model 

The line element for static spherically symmetric spacetimes is given by 

ds 2 = -e v dt 2 + e x dr 2 + r 2 (dd 2 + sin 2 (9d</> 2 ) (1) 

where v{r) and A(r) are arbitrary functions. We study non-radiating rel- 
ativistic spheres with anisotropic stress, and the energy-momentum tensor 
is of the form 

T ab = [iu a u b + ph ab + ir ab (2) 

where fj, is the energy density, p the isotropic pressure, and the projection 
tensor h ab = u a u b + g ab is measured relative to the four- velocity u a . It is 
convenient to express the anisotropic stress in the form 

7T ab = V3S{r) (c a c b - ^h ab ^j 
3 



where the unit spacelike vector c is orthogonal to the fluid four-velocity u 
and IS^r)! is the magnitude of the stress tensor. This representation for 
ir ab is a consequence of the symmetries of the static spherically symmetric 
13j. The quantity S is a useful device to introduce 



spacetimes 



p r =p+2S/V3, p±=p-S/V3 



which are the radial and tangential pressures respectively. Note that for 
isotropic matter S = and p r = p± = p. The magnitude S provides a 
measure of anisotropy. We assume that the fluid four-velocity is comoving. 
This assumption implies that u a — e - "/ 2 ^, c a — e -A / 2 <5°. 

Using and (J2J, the Einstein field equations become 



e 



-A 



— (l-AV-e A ) = /i (3a) 



e 



-A 



— (l-e x + rv') = p r (3b) 



— 2v" - u'X' + + — = p± (3c) 

4 \ r r J 

for static spherically symmetric anisotropic matter. We are using units 
where the speed of light and the coupling constant are unity. The momen- 
tum conservation equation leads to 

{li+p r )v' + 2p' r + ^( Pr - Pl _) = (4) 

for the spacetime This conservation equation is not independent and 
can be generated directly from the field equations 10. We define the mass 
function as 

1 f 

m(r) = — x 2 /i(x)dx (5) 

2 Jo 



following the treatment of Stephani |23|. With the help of (QJ and J5J) we 
can integrate l|3af) . and then get the equivalent system 

\ „ 2m ,„ . 

e~ x = 1 (6a) 

r 

r(r — 2m)v' — p r r 3 + 2m (6b) 
4 

(fl + p r ) V + 2p' r = ~(p r -P±) (6c) 

r 

The system (jSJ has the advantage of being a first order system of differ- 
ential equations, and is linear in the gravitational potential v which sim- 
plifies the integration process. For certain applications it is easier to use 
JBJ rather than the original second order system which is the approach 
that we follow in this paper. We seek explicit solutions to the Einstein field 
equations that describe realistic anisotropic relativistic stars by utilising an 
algorithm that was initially proposed by Maharaj and Maartens |b|. In 
their approach they expressed the field equations as the first order sys- 
tem of differential equations 10 ■ The energy density fi (or equivalently m) 
and the radial pressure p r are chosen on physical grounds. The remaining 
relevant quantities (e^ ,e x ,p±) then follow from the field equations. Note 
that (e^e , p (or m),p r ,p±) are not independent; there are five unknown 
functions and three field equations so that we have the freedom to choose 
any two of the quantities. In this paper we make explicit choices for \x and 

Pr- 
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3 General solution to the field equations 



It is convenient to make the following choice for the energy density 



H = ^ + k + £r 2 (7) 



where j, k and t are constants. The roles of j, k and £ in the physics of the 
model are highlighted in examples considered later. An advantage of this 
form for /i is that it contains particular cases studied previously. Then Ipjl 
yields the following expression for the mass function 

m = i( i+ r 2+ r 4 ) (8) 

with the particular energy density Q. Equation l|6al) gives 

e- A = l-j-^r 2 - e -r 4 (9) 
and the gravitational potential A has been determined. 
With the help of {SJ, we can write (|6bl) as 

l 1 Pr , 



l_j_| r 2_| r 4 r (1 - j - | r 2 _ | r 4> 

rPr j 

1 _ ? _| r2 _| r 4 r(l-j) 

1 3' 5 



(10) 



(l-i)(l-j-|r2-fr4) 
where we have used partial fractions. On integration, (|10|l can be expressed 



as 



7 1 

Ix + ^— lnr+ 7 2 + ln£ (11) 

1 - i i - i 



where In B is a constant of integration and we have set 



h = / : r-r dr 



1 


-j- 


k r 2 t r i 
3' 5' 










fr 


+ | r 3 


1 


-3- 
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h = / — 



At this point we could choose a barotropic equation of state p r — p r (/J-). 
However this is an approach that we intend to follow in future work. In this 
treatment we make a choice for the radial pressure p r which is physically 
reasonable and is a generalisation of earlier studies. We make the choice 

» ^ i?j( 1 -'-r' -£*)('-£)" < 12) 

When j = £ = 0, we obtain the radial pressure postulated by Maharaj and 
Maartens 14 1. For j = 0, we regain the radial pressure of Gokhroo and 
Mehra ^. The form (|12Jl for p r is physically reasonable because p r > 
in the interval (0, R) for relevant choices of the constants, p r = C at the 
centre r = 0, p r — at the boundary r = R, and p r is continuous and well 
behaved in the interval [0, R}. 

The first integral Ji simplifies to 

h = — I (l-iVrdr 



1-3 J V R 2 J 
CR 2 



2(l-i)(n+l)V R 2 ) 
for the choice of l|12f> . To evaluate the second integral I2 we need to consider 
two cases: £ — and £ ^ 0. 
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Case I: £ = 



In this case the integration is straightforward and we obtain 

If k o 



Case II: £^Q 



For this case we let 



and obtain 



h = 



5k 



t (r 2 + 1; 



1 — j + Ml 



5 A' 



I ( r 2 I 5k\ 

36i 5 V ~ 61 ) 



rdr 
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du 



5k 2 



6i' q 1 J + m 



Iniq 



5k ( y/5 



tanh 



10 ^ 21"' Y 1 5 J ' U \qVej \qV5~ 
Hence we can collectively write for both Case I and Case II that 



|ln{l- j - fr 2 } , for £ = 



/ 2 = < 



iln{l-J + |g-f(- 2 + i) 2 } 



2 r + fiC 



(13) 



for I ^ 



The integrals I\ and J2 are given in terms of elementary functions which 
helps in the physical analysis of the model. 



On substituting I\ in Ijllj l we obtain 



(14) 
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for the gravitational potential e v where I2 = Iz{t) has the functional rep- 
resentation given above in (|13[1 for 1 = and I 7^ 0. To match the interior 
solution to the Schwarzschild interior we require that e^ fl) = 1 - 2M/R 
which implies that B = R^J (1 - 2M/R)exp{-I 2 (R)/(j - 1)}. Finally 
the last field equation (|pc|l gives the tangential pressure p± : 



c ( . t 



r 



2 \ n 



t { k o $ a ^ 

H 1-7 r 2 - -r 4 

2 1^3 5 



C 2 A . k 2 * 4 W. r 2 \ 



l-i-^--r» 1 



2(1 -i) 2 V 3 5 J \ R 2 J 



2 / „2 \ 2 ™ 



- 



2nC / fc o t A \* f. r 2 



l-J-^r 2 -^) 1 ' 



2 / o \ »-l 



(1 - i)i? 2 v 3 5 ; v ^ 2 

1 I i , , /?„2 1 (.•,"'2,^4 



+ _^ + fe + ^^- + ir ' + - r 'jj (15) 
where we have used 0, (fT2l and lfH| . The anisotropic factor S'(r) is given 

by 



2 / r 2 \ 2n 



■ ^ 2 ^4 
1 — j — — r r 



2V3 V 3 5 



C /* . k <-) £ 



(1 - i)i? 2 V J 3 5 J \ R 2 J 

+ h{i +k+h2 ){ i+ r 2+ r > )} < w > 

which follows from p2|) and (|15ll . 

Thus we have generated a new class of solutions to the Einstein field equa- 
tions ©■ Collecting the various results given above we can express the 



exact solution as 



= ^ + k + £r 2 (17a) 



C 2 ( . k 



, i 1 - i r* - -r* 1 ^ 

2(1 -j) 2 V 3 5 M R 2 



2 / 2 \ 2 ™ 



- 



(1 - j)R 2 \ J 3 5 / V # 2 / 

+ H^ +k+er2 )(; ]+ r 2+ r 4 )} (17c) 

e" = Br^exp/-^ f i _ \ (17d) 

\W 2(l-j)(n + l)V i?j j 1 ' 

e " = ' r 2_ V ( 17e ) 

J 3 5 

where ^2 (given in jE} ) contains the two cases, I — and ^ 7^ 0. The exact 
solution H17J) represents the interior of an anisotropic star corresponding to 
the energy density /i = j /r 2 + k + lr 2 . Clearly other choices for fi and p r 
will yield new solutions to the field equations; however these choices may 
not correspond to realistic matter or the integrals I\ and I2 may not be ex- 
pressible in closed form. This solution does not have a barotropic equation 
of state p r = p r (/■*)■ To obtain a model with an equation of state we need 
to specify this explicitly when evaluating the integral I\. An equation of 
state is a desirable physical feature which we hope to incorporate in future 
models. Note that the cosmological constant is absent from our model. 
This quantity can be easily included by adding a constant to the energy 
density and the pressure function. 



10 



4 Special cases 



We consider some special cases contained in the new class of solution pre- 
sented in section [3] two of these cases lead to particular models that have 
been studied previously. 

Solution I: j = 

In this case the energy density is given by 

(i = k + £r 2 
and the line element has the form 



1 - -r 2 - 



£ \ _1 

-r 4 J dr 2 + r 2 (d6 2 + sin 2 <9d</> 2 ) (18) 



The particular solution QlSj l was found by Gokhroo and Mehra 8|. Their 
solution is regained when we set k = pa, I = —poK/a 2 . Note that if 
we require fi' < then the constant I < for a monotonically decreasing 
energy density as we approach the boundary r = R from the centre. 

Solution II: j = I = 

For this case the energy density 

p = k 
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is a constant. The line element has the representation 

* - -H-H-*o-^(-;r})* 

1 - K 2 ^j dr 2 + r 2 (d8 2 + sin 2 Odcf 2 ) (19) 

The particular solution l|19|) was found by Maharaj and Maartens 14 1. 
Their solution is regained when we let k = 6M/R 3 . Since fi is constant we 
may interpret this solution as an anisotropic generalisation of the incom- 
pressible Schwarzschild interior sphere; however note that the anisotropy 
factor S(r) ^ everywhere except at the centre r = 0. 

Solution III: k = I = 

In this case the energy density has the form 



M = — 



The line element is given by 
ds 2 = -Bri=J 



X GXP (~2 lD(1 -^- 2(1 - J){ n + 1) i 1 -^ ) d<2 
+ (1 - j)' 1 dr 2 + r 2 (d9 2 + sin 2 dd<p 2 ) (20) 



Even though (|20|l has a very simple form, we believe that it is a new 
anisotropic solution to the Einstein field equations and has not been pub- 
lished before. Since /j, oc r^_ 2 we may relate l|2(Jfl to the results of other treat- 

3 



ments. Dev and Gleiser 



, Herrera and Santos [ll| and Petri 



found 



solutions to the anisotropic Einstein field equations involving fi « r 2 . In 

12 



each of these papers a different set of assumptions to that utilised in this 
paper was used; in our treatment we have chosen a form for the radial 
pressure p r . Therefore their solutions are necessarily different from 12(J[) 
for the corresponding energy density choice \i oc r -2 . 



5 Physical Conditions and Analysis 



One of the original reasons for studying anisotropic matter was to generate 
models that permit redshifts higher than the critical redshift z c of isotropic 
matter 3j. Observational results indicate that certain isolated objects have 
redshifts higher than z c . The surface redshift is given by 



R J 



The critical redshift z c — 2 is the limiting value for the perfect fluid spheres, 
and is attained when 2M/R = 8/9 J. For the range of values falling in the 
interval 8/9 < 2M/R < 1 the redshift is greater than z c ; this phenomenon 
may be explained by allowing for anisotropy. For values of 2M/R close 
to unity, the surface redshift becomes infinitely large. The feasibility of 
higher redshifts for anisotropic matter, in both Newtonian and relativistic 



models, was firmly established by Bondi 



It is interesting to note that 



17| emphasise 



Bondi , Binney and Tremaine |l| , Cuddeford |5j and Michie 
the significance of anisotropies in stellar clusters and galaxies, in addition 
to individual stars. 
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The gravitational potential e A is finite at the centre r = and at the bound- 
ary r — R. The function e A is well behaved in the interior of the relativistic 
star. The gravitational potential e u is continuous and well behaved in the 
interior and finite at the boundary of the star r = R. There is a singularity 
at the centre r = in the potential e v . The singularity in e" is removable 
for a specific choice of parameter values. This singularity is eliminated by 
setting j = which corresponds to the solution of Gokhroo and Mehra 8j . 

The energy density Q chosen describes relativistic stars as we demonstrate 
later. The form of /i oc?'~ 2 (fc = ^ = 0)is usually used in domains where it 
is not possible to use a single equation of state; particularly where the origin 
is excluded, like a body with a constant density core and matter density dis- 



ensity 

1m 



. It is interesting to observe 



tribution around the core going like : 
that the r~ 2 profile in the energy density also arises in isothermal spheres 
in Newtonian configurations that correspond to a Maxwell-Boltzmann gas 
in galactic systems 21 . Densities with j ^ and k ^ are also physi- 
cally reasonable. For example, Misner and Zapolsky |lS| propose that the 
term jr~ 2 models the physical configuration of a relativistic Fermi gas for 
some particular value of the parameter j. Another example is due to Dev 



and Gleiser 



G] 



who suggest that for some particular value of j and k =/= 
the energy density function jr~ 2 + k describes a relativistic Fermi gas core 
immersed in a constant density background. 

The radial pressure p r is continuous and well behaved in the interior of the 
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star. Alsop r > in the interval (0, R), regular at the centre (p r (r = 0) = C), 
and vanishes at the boundary (p r (r = R) =0). The tangential pressure p± 
has a singularity at the centre, but is otherwise well behaved throughout the 
interior of the star and finite at the boundary. The singularity in p± may be 
eliminated by suitable particular choice of parameter values. In general the 
tangential pressure is not zero at the boundary of the star (p±_ (r = R) =^ 0) 
which is different from the radial pressure (p r (r — R) — 0). It is also im- 
portant to observe that the magnitude of the stress tensor 

S = -j=(p r -p±) 

is a nonzero function in general. Hence this class of solutions is gener- 
ally anisotropic and does not have an isotropic limit (the isotropic limit 
results when we set particular values for the constants in our ansatz). It 
is not possible to eliminate S and obtain an isotropic counterpart. This 
means that the model remains anisotropic. An analogous situation rises 
in Einstein-Maxwell solutions modelling charged relativistic stars in which 
the electric field is always present. An example of such a charged star is 



given by Hansraj |9(. 

Figures ^ 12 and|21are illustrations of the behaviour of the energy density 
(U, the radial pressure p r , and the anisotropy factor S respectively, for 
particular chosen values of the constants in the exact solution l|17|) ■ The 
radial distance is over the interval < r < 1 and the boundary of the star 
has been normalised to be r = R = 1. Note that Plot A corresponds to 
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j = case (Solution I), Plot B corresponds to j = i = (Solution II) 
and Plot C corresponds to the general solution l|17fl where j ^ 0, fc ^ 0, 
and £ ^ 0. In Figure ^ for /i, Plots A and B are continuous throughout 
the interval < r < 1; however Plot C indicates unphysical behaviour as 
we approach the centre. This undesirable feature in Plot C arises because 
7^0 and indicates that another solution has to be utilised around the 

onn 

centre in a core-envelope model |6|, llll LUJ • In Figure |U for p r , the radial 
pressure is monotonically decreasing from the centre to the boundary for 
all Plots A, B, and C. In Figure |21 for S, Plots A and B are continuous 
throughout the interval < r < 1; however Plot C indicates a singularity 
as we approach r = 0. We suspect that this singularity in Plot C is related 
to the fact that j ' ^ 0. We observe that the gradient of S is greatest for 
Plot C, corresponding the case for j ^ 0, k ^ 0, and £ ^ 0, as the boundary 
is approached. Hence S(R) has the largest value at the boundary for the 
general solution (|17(l in this case. The behaviour of S outside the centre 
is likely to correspond to physically reasonable anisotropic matter: Plot B 
has a profile similar to the behaviour of the anisotropic boson stars studied 
by Dev and Gleiser p . However the general solution obtained by Dev and 
Gleiser || for the choice /j = jr~ 2 + k is given in terms of hypergeometric 
functions. Our corresponding solution has the advantage of being expressed 
in terms of elementary functions. 
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Figure 1: Energy density /i(r) plots 
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Figure 2: Radial pressure p r (r) plots 
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Figure 3: Anisotropy factor S(r) plots 
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6 Discussion 



It is possible to demonstrate that the solutions found can be utilised to 
discuss the structure of neutron stars and quasi-stellar objects. We write 
the surface density in the particular form /j, s = /io(l — j — t) where (Iq is 
the central density and the constants j, k and I have been chosen so that 
the energy density can be easily expressed in c.g.s. units. Now consider a 
neutron star of radius 10km and surface density of 2 x 10 14 gcm~ 3 . Then 
the parameters (j,q, HqR 2 , 2M/R, surface redshift z = (1 — 2M/R)~ 1 ^ 2 — 1, 
and mass M in terms of the solar masses M P) can be calculated. We 

- - D 

choose values of j and i so that comparison with Gokhroo and Mehra |g] 
is facilitated. The results are given in Table In this category of results 
the surface redshifts range up to 0.566, and masses extend to 2.00Af Q . 
This range of values is consistent with the results of Gokhroo and Mehra 
J|. Hence our solutions yield values for surface redshifts and masses that 
correspond to realistic stellar sources such as Her X-l and Vela X-l. Clearly 
higher values for z and M can be generated by adjusting j and I. 

In this paper we have found a new class of solutions to the Einstein field 
equations for an anisotropic matter distribution utilising the algorithm of 
Maharaj and Maartens 14] . These solutions correspond to the energy 
density fi = jr~ 2 + k + £r 2 and contain particular solutions found previ- 
ously. We note that the term containin g jr ~ 2 is physically important and 



arises in a number of applications 



3 2, m 3 



Our results indicate that 
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3 


I 


Mo x 10 14 


[i R 2 


2M/R 


z 


M(M Q ) 








1.48 


.015 


.124 


.068 


.42 


.001 


.1 


1.65 


.017 


.130 


.072 


.44 


.002 


.2 


1.86 


.019 


.136 


.076 


.46 


.003 


.3 


2.13 


.021 


.145 


.081 


.49 


.004 


.4 


2.49 


.025 


.156 


.089 


.53 


.005 


.5 


3.00 


.030 


.172 


.099 


.58 


.006 


.6 


3.77 


.038 


.196 


.116 


.66 


.007 


.7 


5.07 


.051 


.237 


.145 


.80 


.008 


.8 


7.73 


.077 


.321 


.214 


1.09 


.009 


.9 


16.32 


.163 


.592 


.566 


2.00 



Tabic 1: Densities and redshifts for neutron stars 



anisotropic solutions for the physically reasonable energy density /j, oc r~ 2 
can be generated with the simple solution generating mechanism of Ma- 
haraj and Maartens 14 1 . Our ongoing investigations indicate that a general 
class of anisotropic models are possible, for different choices of fi, such that 
the desired limit fj, oc r~ 2 is regained as we approach the boundary. This 
work is in preparation. Observe from (|16|) that the anisotropy factor S is 
nonzero in general in the interior of the star. This means that the exact 
solution l|17|) remains anisotropic and does not have an isotropic limit. We 
would need to use another approach of integrating the anisotropic Einstein 
field equations than the algorithm used in this paper, if an isotropic limit 
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is to be contained in the stellar model. 
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